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Stark-induced electric dipole amplitudes between states of the same nominal parity can be impor- 
tant in experiments to observe parity nonconservation in atoms. The Stark-induced El amplitudes 
are expressed in terms of an irreducible spherical-tensor decomposition. This formalism is applied to 
the specific case of transitions between hyperfine sublevels of a single atomic state. It is shown that 
in the ground states of alkali atoms, such transitions are suppressed by many orders of magnitude 
relative to naive expectations. 
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I. INTRODUCTION 

o 

The first definitive observation of a nuclear spin-dependent (NSD) contribution to atomic parity nonconservation 
C/3 ■ (PNC) was recently reported JjJ. This effect arises primarily because of an electromagnetic interaction between 
atomic electrons and the nucleus. In particular, because of PNC interactions within the nucleus, the nucleus acquires 
a P-odd electromagnetic multipole moment known as the anapole moment . Measurements of atomic NSD-PNC 



amplitudes can thus be interpreted in terms of the nuclear anapole moment. The nuclear anapole moment itself can 
be related to the strengths of various hadronic PNC couplings. These couplings are in general poorly determined, and 
the available data are in some cases contradictory. Thus, there is considerable interest in additional measurements of 
atomic NSD-PNC, as a means to probe hadronic neutral-weak interactions 

NSD-PNC effects are extremely small, and their observation requires both extreme statistical sensitivity and careful 
control over minute systematic effects. Gorshkov et al. have proposed a novel technique for observing NSD-PNC [[|. 
. They suggest using a strong microwave electric field to drive PNC-induced electric dipole (El) transitions between 
hyperfine sublevels of a single atomic state. Interference between this El transition amplitude and a (parity-allowed) 
y— ( , Ml amplitude driven by a weak microwave magnetic field (coherent with the electric field) gives rise to a PNC 
observable. This technique offers the promise for unprecedented statistical sensitivity. However, Gorshkov et al. give 
only a limited and qualitative analysis of possible systematic effects in the measurement they propose. 

Further consideration of the proposal of Ref. || has led us to investigate the phenomenon of Stark-induced (SI) 
r/3 | El transitions between hyperfine sublevels. These arise in the presence of a DC electric field, which mixes states of 
opposite parity into the initial and final states; thus El transitions can be induced between one nominal state and 
the admixed component of the other state . Such transitions can be of importance in PNC measurements for two 
reasons. First, the presence of stray electric fields leads to SI-E1 transition amplitudes between hyperfine sublevels, 
which can mimic the effect of the PNC-induced El amplitudes; the size of the SI amplitudes relative to the PNC- 
induced amplitudes determines the extent to which stray fields must be controlled. Second, by deliberately applying 
a large electric field, the SI-E1 amplitude could in principle be used as the primary parity-allowed amplitude, against 
which the PNC amplitude can interfere |t|]. This technique gives excellent control over a variety of systematic effects, 
i and indeed is the principle of the measurement of Ref. ||J], in which the optical 6s- 7s transition in Cs was studied. 

We show here that in the particular case of alkali atoms (as discussed in Ref. ||), the SI-E1 amplitudes between 
hyperfine sublevels of the ground state are suppressed by many orders of magnitude relative to naive expectations. This 
conclusion means that it is impossible in practice to use the SI amplitudes as the primary parity-allowed amplitude 
in experiments of the type proposed by Gorshkov et al. On the other hand, it also means that systematic effects due 
to stray electric fields in such experiments are negligibly small under reasonable conditions. 

II. GENERAL ARGUMENT FOR THE SUPPRESSION OF SI-E1 HYPERFINE TRANSITIONS 

The SI-E1 transition a — > b in a constant field E and time-varying field e has amplitude 
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W ba = J2(- 1 ) 9+9 ' A i,i'E- q e- q ', (1) 



where 



A^ = ^ r ^ r + -- K ^ r ) , (2) 

and the subscripts q, q' refer to spherical vector components. 

The SI-E1 transition can be usefully understood as a special case of a two-photon transition, with one photon 
at zero frequency. For electric dipole interactions, each photon has orbital angular momentum L = with respect 
to the atom. Thus, since each photon has spin S = 1, the two-photon system must have total angular momentum 
Jtot — <Stot — 0,1, or 2. These three cases correspond to scalar (ast), vector (/3gt) , and tensor (7st) transition 
amplitudes, respectively. We are primarily interested here in transitions between sublevels of the s 1 / 2 ground state 
of an alkali atom. For this case, each of the three possible amplitudes is strongly suppressed (or vanishes) relative to 
the naive estimate of the amplitude: A ~ e 2 a^/TZ ~ Oq (where do is the Bohr radius and 1Z is the Rydberg constant). 
These suppressions can be understood in terms of selection rules for two-photon transitions, which were derived in 
Ref. We begin by considering these (approximate) selection rules for the specific case of hyperfine transitions, and 
estimate the residual values of the SI-E1 amplitudes. 

First, we note that the scalar amplitude ast vanishes for hyperfine transitions. This is a trivial consequence of 
the fact that the scalar amplitude cannot change any angular quantum numbers, and thus cannot contribute to any 
transitions between sublevels of the same state. This argument fails when the principal quantum number changes; 
thus, for the well-studied 6s- 7s transition in Cs, ast (6s — 7s) = 269 ||. (This is somewhat larger than the 
naive estimate given above because the radial matrix elements in Cs and other alkalis are ~ 5ao, and the energy 
denominators are ~ 0.172..) 

Next, we show that for hyperfine transitions, the vector amplitude fist is suppressed by a factor ~ — Za 2 , where 
a is the fine-structure constant. This suppression can be understood as a consequence of Bose statistics for photons. 
For Stot=l, the spin component of a two-photon wavefunction is antisymmetric under interchange of the two photons. 
Thus, the space-time component of the wavefunction must also be antisymmetric under interchange, in order to 
satisfy Bose statistics. In the electric-dipole approximation, where the spatial part of the wavefunction is constant, 
the total wavefunction must therefore vanish if the photons have the same energy. (This statement is closely related 
to the Landau- Yang theorem of high-energy physics, which states that the decay of a J=l particle into two photons 
is forbidden |lC],[Il|]; this result has also been derived for atomic two-photon transitions B.) In the limit where the 
frequencies of the two photons are close but not identical, the vector part of the two-photon transition amplitude 
does not vanish, but is suppressed by a factor ~ ^I^" 2 > where is the energy of the ith photon, and E v — E r is 
the difference in energy between a virtual intermediate state and a real intermediate state |ll|. For the case of SI-E1 
transitions between hyperfine sublevels of the same state, Tioj\ — huj2 = hu> = Bhfe (the hyperfine splitting energy), 
E v — Ehf s (or 0), and E r = E c \ (the energy splitting between electronic states). Thus, for the case of interest here, 
/3st is suppressed by an additional factor ~ ~ ^Za 2 (see e.g. fL3fl). Clearly, this suppression is absent for the 
vector component of optical SI-E1 transitions such as the 6s- 7s transition in Cs. 

For the particular case of transitions between hyperfine sublevels of the ground state of an alkali atom, the vector 
amplitude is suppressed by an additional factor of ~ Z 2 a 2 . In the absence of both electron and nuclear spin, /?st 
must vanish for any s — s' (0 — 0') transition. However, the addition of electron spin alone is sufficient to allow this 
amplitude, since then the transition is of the type h ~ h ■ (Nuclear spin alone will also induce a nonzero value of 
(3st, but this is a much smaller effect.) Since for any s%/2 — s' x , 2 transition /3st explicitly relies on the presence of 

electron spin, it must be that /3g t is suppressed relative to the naive estimate by ~ e E -e ~ 1^7 ~ Z 2 a 2 , where Ef s 
is a fine-structure energy splitting (see e.g. (llj). This suppression is indeed present for the 6s- 7s transition in Cs [0, 
where the ratio of vector to scalar amplitudes is @st/ a St = 0.1 [H. For the hyperfine transitions of interest here, the 
vector amplitude is thus suppressed relative to its naive value by an overall factor ~ ^Z 3 a . 

Finally, we consider the tensor part 7st of the SI-E1 amplitude. For s — s' transitions this amplitude vanishes even 
in the presence of electron spin (i.e., for & \ — \ transition); its presence relies explicitly on the presence of nuclear 
spin. Thus, 7st must be suppressed relative to the naive value by a factor ~ e E -e ~ "bT ~ m^^® 2 - us ma ^ e 
this statement more explicit by noting that, for any nuclear spin / > 1, nuclear spin alone is sufficient to induce a 
non-zero tensor amplitude. (For all stable alkali-atom nuclei, / > 1.) Since electron spin is neither necessary nor 
sufficient to produce a tensor amplitude, we expect that its presence should not significantly alter the magnitude of 
7St- Rather, the tensor amplitude can be understood from the hypothetical case of an atom with nuclear spin, but 
no electron spin. In this case, the only relevant hyperfine splitting is that due to the interaction of p-states with the 
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nuclear magnetic dipole and electric quadrupole moments. Thus, 7st is suppressed relative to its naive value by a 
factor of ~ Eh te( p > , This argument holds even for optical s — s' transitions such as the 6s-7s transition in Cs; for this 
reason, the amplitude 7st is negligibly small for this transition, and is always ignored. We note in passing that the 
suppression of 7§ t is not present at all for transitions between hyperfine sublevels of any electronic state with J > 1. 
However, we confine ourselves to the particular case of hyperfine transitions in alkali-atom ground states. 



III. EXPRESSIONS FOR SCALAR, VECTOR AND TENSOR POLARIZABILITIES 



We now discuss specific calculations of the SI-E1 amplitudes. For our purposes it is convenient to rewrite the 
general second-rank tensors A q . q i and E- q E- q ' of ( II) and (|J) in terms of their irreducible tensor components Aq and 
(E ® £)-q- We use the standard transformations 13| , e.g. 

A« = (-l)<5x/2ifTI£ ( \ _\, K Q ) A M ,, (3) 

q,q> 

and 

so that the transition probability ([l]) is written as the contraction of irreducible tensor components (as was done for 
the general case of two-photon transitions in Ref. || ) : 

W ba = J2(-V QA Q( E ® £ ^Q- ( 5 ) 

K,Q 

We further rewrite the irreducible tensor component Aq in terms of A K , the reduced amplitude of rank K : 

aK _ (_i\Jb-M b ( J b K J a \ K . . 

A ® - ( 1} ^ -M b Q M a ) [b} 

In the general case the transition amplitude is fully described by three independent reduced amplitudes with K = 
0,1,2. 

The SI-E1 transition amplitude has traditionally been written in terms of the quantities ast, /?St) and 7st, such 
that 

W ba = a st F, ■ e + p st • E x e + ^ t k (^E t e k + ^E k e z - ^E • eS hk ). (7) 
The quantities ast, Pst, and jst are simply related to the reduced tensor amplitudes in Eq. (^); for instance, 

From here on we will use the reduced amplitudes A K rather than ast, /?st, and 7st, since the irreducible tensor 
formalism makes it possible to write expressions in a general form for all K . 

At this stage let us assume that there is no nuclear spin and levels are described by quantum numbers J;, Mj. Then 
we can replace the matrix elements of r q in (Q) in terms of reduced matrix elements, e.g. 

(nK|a) = (-l)^(n|H|a)(_^ J £ ) . (9) 

Using Eqs. (|2|), (^), (^J), and (^|), we obtain an expression for A K , with the sum (over q, q', and M n ) of the product 
of the three 3j-symbols reduced to a 6j-symbol: 



A- = V2Z+lE(-l)^{^ f j\ 



V\\r\\!;u^\\r\\a)[ T J—+ [ i^—). (10) 
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A. SI-E1 amplitude for nonzero nuclear spin 



We now include quantum numbers F, M and I, First we can simply substitute J, M — > F, M in (|Io|): 
A K = V2iTTI$>l)^ { 5 F > £ } x 

(n,F n \\r\\b,F b ){n,F n \\r\\a,F a ) (—1— + ^!—) . (11) 

\ £j b — &n &a — &n J 

The dependence of the reduced matrix element on the quantum numbers F is described by the expression: 

(n,F n ||r||a,F ) = {-l) J ~ +I+F « +1 ^ {2F n + l)(2F a + 1) j £ J [ } (n||r||a), (12) 



so that 



A K = ( _ l)2 Z+2F a +F 6 ^ (2Jf + 1} (2Fn + 1)(2Fq + 1} ^(_l)2J„+F„ (2Fn + 1} 

n 

F R 1 H f F f] f / F fl / 1 (— lV 

[n\\r\\b)(n\\r\\a) 



K 1 F a j \ F b J b 1 j \ F a J a 1 J v 1 \E b -E n E a -E, 

Up to this point, our expressions have been perfectly general. From now on, we introduce expressions that are specific 
to single valence-electron atoms. In particular, we write the reduced matrix elements of r explicitly in terms of the 
single-electron angular momenta / (orbital) and j (total) : 



(n||r||a) = (-iy- +l ^W{2j n + l)(2j a + l)l an { f f 

^ Ja l a J- 



Rn,a, (13) 



where l an = m&x(l a ,l n ) and R n , a is the radial integral. Combining these two equations gives 



A K = -V(2K + l)(2F a + l)(2F b + l)(2j a + \){2j b + 1) £ C$ Mn D* , (14) 

n 

Cj£ ndntln = (-l) F "+ F >+^(2F n + l)(2 Jn + 1) x 

' On \ \ { In jn \ \ /ir\ 
la 1 j\ib /fe 1 J ' 1 j 

(16) 



/ F n 


F b 




/ jn 


F n 


I 




1 


A} 


\F a 


ja 


1 



Dn = (—l) lan+lb " \flanhnRn,aRn.b ( ~= = h 



(-1) 



Note that £)„ oc E a — E b ; this is the suppression discussed earlier, due to Bose statistics for photons. 

B. Summation over F n and j n 

The general arguments for the suppression of the hyperfme SI-E1 amplitudes suggest that there are cancellations 
in the sums over quantum numbers F n and j n . To see this explicitly, note that in the nonrelativistic limit, does 
not depend on j n or F n . Application of the Racah-Eliot relations fig] to the sums over F n and j n gives: 



Cp n , jn ,i n = (-l) F « +I+K+ ^(2 Jn + 1) x 

If; „• 11 r; „• I 1 

(17) 



K ja jb \ / K ja jb \ / In jn | \ / In jn \ 
I F b F a j \ j n 1 1 j\j a l a 1 J \ j b l b 1 / ' 

E <**.«. = (-ir™*-**-* f £ £ } { f„ ! r * } . <«> 

Note that both expressions turn to zero for K > j a +j b - this is equivalent to the vanishing of a | — ^ tensor amplitude 
in the absence of nuclear spin. In addition expression (^) turns to zero for X > l a + l b . This is equivalent to the 
vanishing of s — s' vector (tensor) amplitudes in the absence of electron (nuclear) spin. 
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Equations ( |l7| ) and ( |l8| ) make it possible to simplify the general expression ( |l4| ) . If we neglect the dependence of 
D% on F n , we can use the sum (|l7|). If we also neglect the dependence of D% on j„, we can use the sum (Q. This 
yields the following expressions: 

A" = -V{2'< + + 1)(2F, + l)(2j„ + l)(2j„ + 1) | * j? * 

d-^^w*. + •> { I \ f } { 1 1 } { i t f \ D - • < w > 





3a 


h) 


[ K l a lb \ 




F b 


F a j 


1 5 jfc ja J 



= - V(2X + l)(2F a + 1)(2F 6 + l)(2j a + l)(2j 6 4 

^ ( _ 1) F.+/+i. + A-W.-i/ 2 j * *6 W (20) 

where is given by (|l6|). The first of these equations is correct up to the hyperfine corrections, while the second is 
correct only up to the fine structure corrections. 



IV. SI-E1 AMPLITUDES FOR HYPERFINE TRANSITIONS IN ALKALI ATOMS 

We now apply these formulae to the F a = I — i — > Fj, = I + h transitions of the ground states noSi/2 of alkali 
atoms. As we have seen above, for such transition both vector and tensor amplitudes are strongly suppressed. Here 
we want to obtain accurate estimates of these amplitudes. 

It is known that for all alkalis 

|-Rn p,n s| ^ \Rnp,n s\, Tl > Uq, (21) 

Thus, the dominant contribution to the sums (|l4|), (|l9| ) and (^o|) comes from the p shell with the principle quantum 
number no. 



A. Vector amplitude 

For the case of interest where l a = h = 0, the nonrelativistic expression (^0|) gives zero value for the amplitude 
A 1 ; to get a correct value, we must use equation (|l9|) which includes fine-structure corrections. In order to avoid 
cancellation in the sum over j n we have to take into account the dependence of the factor £)„ on j n , which can be 
written as: 

Di = D l n + (jn - InK, (22) 

where Z)* and d\ do not depend on j n . We have seen that the contributions of j n = i and j n = | cancel each other 
for I)* . Then obviously for d x n they should double each other, and we find: 

^^(-ifvWnjr? I+ J }}£C (23) 
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lip 



Note that the sum runs over states np and does not include summation over j n . 

Now we have to find d\ p . It follows from ( p2| ) that d\ p = a — a and it is not zero because of the spin-orbit 
interaction H so . There are diagonal and off-diagonal contributions of H so to d^ p : 

• The fine-structure splitting changes the energy denominators in (|l6|). If we define A SOl np — Fnp^^ — ^ n Pi/2 ^^d 
note that for the hyperfine transition E — E a = A noS (I + 5), where A naS is the hyperfine constant of the ground 
state, we get for d\: 



-(2J+l)i£ P|no S ( ; noS _ S °™ . (24) 
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The spin-orbit interaction changes the radial integrals in (|l6|): 



c p (np\H so \n'p)R n , p , noS 

ott, npj ,n s — p ■ y zo ) 

^np ^np' 



Parametrically both mechanisms give the same smallness for , but numerically for alkali atoms the off-diagonal 
correction is suppressed. Indeed, for the dominant term n — hq, the correction to the radial integral is numerically 
small because of the relation (|2l|). 

Thus the vector amplitude can be written: 



A A 

^iios^so.np (OR) 
n p,n s/p _ rji \ 3 - \ zo ) 



B. Tensor amplitude 



Here both expressions (|T^) and ( pp| ) turn to zero and we have to return to the expressions (|l4|)-(|l6|). Because of the 
relation ( [21] ) we again can restrict the summation to the shell nop. Then the sum runs over F — I ± A and j = A , | . 
According to (|l7| ) the first sum vanishes because j a = jb = \, while (|l8| ) shows that the second sum vanishes because 
l a = l b = 0. That means that all four coefficients Cj^ * differ only by their signs. So, it is sufficient to calculate, 
for example, Cz L i : 

I Q , Q ,U 



1 1 

2 ' 2 ' 



:(_i)2/+i 7 



i-h i 



l 



a i-i i 

2 J 2 
J- i i 1 



^-1 ' 

A A 1 

2 2 



(27) 



It is clear that in order to avoid cancellations we have to include hyperfine corrections to . Again, as in the case 
of vector amplitude, we can neglect the correction to the radial integrals and consider only corrections to the energy 
denominators. For the «oP3/2 state the quadrupole hyperfine structure should be included. That gives the following 
result: 



A 1 



1)^5/(7+1) 



'S[ 

2 A 



R 



2 



1 

I - 



{E„ 



E Y 



-A 



I - 



2 1 



"0P3/2 ^ (1(21-1)) 



(28) 



where A naPj and B noP3/2 are the hyperfine constants for the levels noPj ■ 



C. Numerical results 

Equations ( p6| ) and ( p8| ) show that, in agreement with our general arguments, the vector amplitude A 1 ~ 
A h fa,, 1Qa a ^."op q 3 x ^Z 3 a 4 , while the tensor amplitude A 2 ~ Ah ^-" 0P ag cx ^Za 2 . This means that for light 
alkalis the tensor amplitude dominates; however, the vector amplitude grows more rapidly with Z and for heavy 
alkalis may even be larger than the tensor amplitude, since the hyperfine constant of the ground s state is much larger 
than the hyperfine constants of the p states. 

To obtain numerical values of these amplitudes we need to know the radial integrals R noP ,n s and hyperfine constants 
A naS and A noPj . They are given in Table | together with our results for the amplitudes A 1 and A 2 . The uncertainty 
in these values arises primarily from the two approximations we have made, namely: neglecting intermediate p-states 
with n > n in the sum of equation (|l4j), and neglecting the j-dependence of radial matrix elements. We estimate 
that these approximations introduce errors of < 10%. 

As stated in the introduction, the motivation for this calculation was the effect of the SI-E1 amplitudes in ex- 
periments to measure PNC in hyperfine transitions. In order to make the Si-El amplitudes comparable to the 
parity-allowed Ml transitions between hyperfine sublevels (and thus to use the SI-E1 amplitudes as the primary 
parity-allowed amplitude), it is necessary that A ~ \±b- Even for Fr, where A attains its largest value, this is achieved 
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only for E ~ 5 x 10 9 V/cm! Such a prospect is clearly unrealistic. On the other hand, it is also instructive to 
compare the magnitude of the SI-E1 amplitudes to the PNC-induced El amplitudes. Using the values for the PNC 
amplitudes for hyperfine transitions from Ref. || (kDc s ~ —3 x 10 _13 eao, kDk ~ — 1 x 10~ 14 eao), we find that in 
both cases the largest Stark-induced amplitude is comparable to the PNC-induced amplitude only when the electric 
field E ~ 3 V/cm. Since stray fields well below this level are easily controlled, systematic effects due to uncontrolled 
SI-E1 amplitudes in experiments of the type proposed by Gorshkov et al. should not be a serious problem, even for 
the lighter alkalis. 

We wish to thank D. Budker and S. Porsev for helpful discussions and useful comments. 
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TABLE I. Input data and numerical values for SI-E1 amplitudes. Unless otherwise noted, data 
for A sp , A so , noP , A nQ3 , A noPl/2 , B noV3/2 , and R nQ s, noP are taken from Ref. jl(J. 







7 Li 


23 Na 


39 K 


85 Rb 


133 Cs 


221 pr 


I 




3 


3 


3 


5 


7 


5 a 




2 


2 


2 


2 


2 


2 


^sp 


(cm" 1 ) 


14904 


16965 


13014 


12698 


11456 


13081 b 




(cm- 1 ) 


0.34 


17.2 


57.9 


238 


554 


1687 b 


s 


(GHz) 


0.402 


0.886 


0.231 


1.01 


2.30 


6.21 a 




(MHz) 


45.9 


94.3 


27.8 


120.7 


292 


811 b 




(MHz) 


-3.06 


18.7 


6.1 


25.0 


50.3 


66.5 b 


Bn p 3/2 


(MHz) 


-0.2 


2.9 


2.8 


26.0 


-0.4 


-260 b 


RriQ s,tiqp 


(o ) 


4.05 


4.29 


5.06 


5.03 


5.50 


5.11 c 


A 1 


(10- 6 ag) 


0.0085 


0.72 


2.0 


70 


940 


2900 


A 2 


(10- 6 ag) 


20 


27 


20 


185 


1050 


1350 



a Ref. Hi 
b Ref. |J 
c Ref. ||f 
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